The Strange Quark Contribution to the Proton's Magnetic Moment 
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We report a new determination of the strange quark contribution to the proton's magnetic form 
factor at a four-momentum transfer = 0.1(GeV/c)^ from parity-violating e-p elastic scattering. 
The result uses a revised analysis of data from the SAMPLE experiment which was carried out at the 
MIT-Bates Laboratory. The data are combined with a calculation of the proton's axial form factor 
Ga to determine the strange form factor G'l/((5^=0.1)=0.37±0.20±0.26±0.07. The extrapolation 
of Gil to its = limit and comparison with calculations is also discussed. 
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In 1988, Kaplan and Manohar jl| showed that unique 
information about the contribution of sea quarks to 
ground state nucleon properties, such as its spin, charge 
and magnetic moment, can be gained through study of 
neutral weak probes of the nucleon such as neutrino- 
nucleon scattering. Soon after, McKeown and Beck Q 
outlined a possible program of experiments in parity- 
violating electron scattering that, when combined with 
existing measurements of nucleon electromagnetic form 
factors, would allow the identification of possible strange 
quark contributions to the proton's charge and mag- 
netism. Experiments have since been carried out H, |y, 
0,0 and proposed 0,^3 several electron scattering fa- 
cilities, and many theoretical calculations of the strange 
quark components of these neutral weak matrix elements 
have appeared in the literature. 

In this paper we present a new determination of the 
strange quark contribution to the proton's magnetic form 
factor, G%j, from the parity- violating e-p elastic scatter- 
ing data of the SAMPLE experiment, which was carried 
out at the MIT-Bates Laboratory in 1998. Preliminary 
results were presented in in the intervening period 
two measurements with a deuterium target were per- 
formed along with a more detailed analysis of the hy- 
drogen data. The deuterium data are relatively insensi- 
tive to the strange vector matrix elements, but provide 
the first experimental information about the nucleon's 
neutral weak axial form factor as seen by an electron 
probe, G^, which is necessary for a reliable extraction 



of G%,j. Preliminary results using the first of the two 
deuterium measurements indicated a relatively small 
contribution from strange quarks to the proton's mag- 
netic moment but an unexpectedly large deviation of the 
axial form factor from theoretical expectation as com- 
puted in 0. These data have been further analyzed, 
and a new experiment was carried out at lower momen- 
tum transfer, and both measurements are now in good 
agreement with the calculation. The deuterium data are 
discussed in a separate report [l^ . Here we re-evaluate 
the implications of the hydrogen data in light of this new 
conclusion, using the calculation of for the axial form 
factor. We also assess its implication for the contribution 
of strange quarks to the proton's magnetic moment using 
a theoretical extrapolation of to its static limit. 

In parity-violating elastic e-p scattering, the asymme- 
try in the scattering cross section with respect to the 
incident electron's helicity is 
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e' = a/t (1 -I- t) (1 — e"^). At the backward-angle kine- 
matics of the SAMPLE experiment, the e-p asymmetry is 
dominated by the contribution from the magnetic neutral 
weak form factor, Gfj. When combined with measure- 
ments of the proton and neutron electromagnetic form 
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factors, the strange quark component can be extracted 
through the relation 

GIj = (1 - 4sin2 Ow) (1 + B^y) GIj~{1 + R^) GIj-GI . 

_ (2) 

The radiative corrections |l3j represent (small) con- 
tributions from higher order processes. While the SAM- 
PLE measurement is dominated by this term, it is also 
sensitive to the proton's neutral weak axial form factor 
G\. At tree level, isospin symmetry relates G\ to the 
axial vector coupling that enters neutron decay. Small 
corrections of order (—10)% are generated by strange 
quarks Il4'|. Electroweak radiative corrections introduce 
a more substantial, 0(— 50%) effect. While electroweak 
corrections that renormalize the individual quark axial 
vector currents can be reliably computed, other effects 
that involve strong and weak qq correlations, such as 
Z— 7 box graphs and the nucleon "anapole moment" [isj . 
present a theoretical challenge. Predictions for these 
corrections, including estimates of the theoretical uncer- 
tainties obtained with chiral perturbation theory (xPT), 
were obtained in j3| and updated in The results 

of the SAMPLE deuterium measurements are consistent 
with these predictions, and in what follows, we use the 
value for G\ obtained in to interpret the results of 
the SAMPLE e-p measurement. 

The SAMPLE e-p measurement was carried out at the 
MIT-Bates Linear Accelerator Center in 1998. A beam 
of 200 MeV circularly polarized electrons was incident on 
a 40 cm liquid hydrogen target, and Cerenkov light from 
backward-scattered electrons were detected in an array 
of ten photomultiplier tubes after reflection from ellip- 
soidal mirrors. The yield in each photomultiplier tube 
was integrated over the 25 /usee long beam pulse and 
sorted by beam helicity state. The measured asymme- 
try was computed as the difference between the yields 
in the two helicity states over the sum, with corrections 
coming from helicity correlations in the beam, and from 
dilution factors associated with the beam polarization 
(36.2±0.1%), electromagnetic radiative corrections, and 
electromagnetic background (~30% of the yield). Pre- 
liminary results were published in , in which additional 
details of the experimental method can be found. 

Subsequent refinement of the data analysis, along with 
development of a GEANT-based Monte Carlo simula- 
tion [l7| of the full experimental geometry, has revealed 
three corrections, all of which act to increase the magni- 
tude of the experimental asymmetry. First, the electro- 
magnetic radiative corrections were recomputed within 
the context of the simulation, whereas in jSj they were 
computed at the central kinematics of each detector, re- 
sulting in a 4% increase in the dilution factor. In both 
cases a spin-dependent modification to the formalism 
of Mo and Tsai was used to compute the radiative 
effects. In the simulation, scattered electron events were 
generated uniformly in energy, angle and along the length 



of the 40 cm target. Energy loss due to ionization and 
collisions in the aluminum entrance window to the tar- 
get, and in the thickness of liquid hydrogen upstream of 
the randomly chosen interaction point, was accounted for 
before computation of the scattered electron kinematics. 
Each scattered electron was assigned a cross section and 
a parity-violating asymmetry, and propagated through 
the target exit windows and the scattering chamber. A 
detection efficiency based on the velocity of the outgo- 
ing electron and the path length of the event's track in 
the Cerenkov medium was combined with the computed 
cross section as an event weight. The radiative correction 
factor of approximately 1.13 was evaluated separately for 
each detector module, and was computed as the ratio of 
the (weighted) asymmetry without and with the radia- 
tive effects included. 

Secondly, a background associated with threshold 
photo- pion production, which had been neglected in j^, 
was evaluated using the GEANT simulation. Such pro- 
cesses contribute to the detector yield through their de- 
cay products, but have a negligible parity- violating asym- 
metry |23|. The 7r° (tt"*") channel was modeled based 
on data from ([13 )■ The 7r+ production yield was 
found to be consistent with experimental observation of 
an exponential tail in the detector signal corresponding 
to the arrival of decay products of secondary muons. The 
net additional dilution factor coming from the pion back- 
ground was 1.04. 

The third modification to the previous analysis was in 
the treatment of the background coming from charged 
particles that were not blocked by shutters placed in 
front of the photomultipher tubes. As discussed in |^, 
the net measured "shutter closed" background asym- 
metry was consistent with zero, but the detector-by- 
detector distribution appeared to have a nonstatistical 
component, and a systematic error accounting for the 
nonstatistical behavior was added to the experimental 
uncertainty for each individual detector before combin- 
ing them. Subsequent analysis revealed that this shut- 
ter closed distribution had a (/)-dependence which fit the 
function /(0) — Aq+Ai cos (2(/) -I- 0o) with a significantly 
better x^/d.o.f (8.0/7) than the presumed flat distribu- 
tion (33.3/9). The OPEN shutter data did not show 
such behavior with statistical significance. The co- 
efficient, — 0.06±0.71 ppm after all dilution corrections, 
was then subtracted from the OPEN asymmetry and its 
uncertainty added in quadrature. This method produced 
a 5% larger result than the method used in 5]. The three 
effects combined result in an experimental e-p parity- 
violating asymmetry of ^23| 

A(Q2 = 0.1) = -5.61 ± 0.67 ± 0.88 ppm . (3) 

The Monte Carlo simulation was also used to determine 
the appropriate theoretical asymmetry to which the data 
should be compared. Averaging over detector and target 
length acceptance effects results in an approximately 3% 
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smaller theoretical asymmetry, 

A(Q2 ^ 0.1) = -5.56 + 3.S7Gli + 1.54G^(^=^) ppm. 

(4) 

The small isoscalar component of has been absorbed 
into the first term. In the model, dipole form factors were 
used for and G"^ , and the Galster parameteriza- 

tion HI for'G|. 

The two SAMPLE deuterium measurements were also 
analyzed using the GEANT simulation, as described 
in |l2|. Both measurements are in agreement with the 
theoretical prediction for the asymmetry using the elec- 
troweak radiative corrections of |lll| in the computation 
of G^'^=^)(Q2). The results for the two 200 MeV data 
sets, along with the computation of are shown in 
Figure Has 1-cr bands in the space of G%j vs. G\^^~^\ 
Both the overlap of the hydrogen data and the calcu- 
lation, and the overlap of the two data sets, are shown 
as ellipses demonstrating the good agreement be- 

tween the deuterium data and the theoretical expectation 
of G^('^=^' = -0.83 ± 0.26. Using this value results in 

Gm(Q^ = 0.1) = 0.37 ± 0.20 ± 0.26 ± 0.07 (5) 

where the uncertainties are statistical, experimental sys- 
tematic and the uncertainty due to electroweak radiative 
corrections, respectively. Inclusion of the SAMPLE deu- 
terium data in the extraction of G|j produces essentially 
the same result. Our new result is slightly shifted in the 
positive direction relative to, but consistent with, the 
analysis in Ref. W\. 

The SAMPLE results for G^j, together with the con- 
siderable body of theoretical work on this topic under- 
taken over the past decade, has substantially sharpened 
our picture of strange quark contributions to the nu- 
cleon's electromagnetic structure. Here, we summarize 
the insights provided by the various theoretical studies of 
strange magnetism (for a more extensive review, see j25jV 

The most direct, first principles approach in QCD re- 
lies on lattice simulations. Such calculations are partic- 
ularly challenging for operators such as S7pS that do not 
connect with nucleon sources, since one must sum dis- 
connected quark loops over all lattice sites. Although 
a complete calculation is prohibitively time consuming 
and expensive, numerically tractable estimates can be ob- 
tained by approximating the sum with so-called Z2 noise 
methods [23. A few calculations have been carried out, 
but with little consensus on the value of G\j{Q'^). The 
calculation in jl^l that relied on a relatively small sample 
of gauge field configurations and large number of Z2 noise 
vectors yielded a statistically significant, non-zero, neg- 
ative value for G\.j over a range of momentum transfer, 
with /is — G|j(0) = —0.36 ± 0.20. In contrast, the com- 
putation reported in |28l | , which drew upon a much larger 
gauge configuration sample and smaller sample of noise 
vectors, found = 0.05 ± 0.06. A lattice-based com- 
putation by the Adelaide group is similarly consistent 
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FIG. 1: Uncertainty bands of Gl, vs. 6^*^=^' at Q^^O.l 
(GeV/c)^ resulting from the two 200 MeV data sets of the 
SAMPLE experiment. Also shown is the uncertainty band 
of the theoretical expectation G^^'"^^^' as computed by pl| 
extrapolated to the same momentum transfer. The smaller 
ellipse corresponds to the 1 o overlap of the hydrogen data 
and the theoretical prediction, the larger one the l-cr overlap 
of the two data sets. 



with zero [29j, and our experimental result tends to fa- 
vor these more recent calculations. One anticipates that 
further refinements in lattice methods, such as the use 
of chiral fermions and partially-quenched or unquenched 
computations, will eventually lead to a consensus for the 
precise value of G\[ . 

Alternative treatments have been pursued using dis- 
persion relations and chiral perturbation theory (xPT). 
In principle, these methods are model-independent since 
they rely on general properties of QCD, such as analyt- 
icity and causality (dispersion relations) or the approx- 
imate chiral symmetry of QCD (^PT) to organize in a 
hadronic basis various contributions to G|j. In either 
case, one must rely upon experimental input to determine 
these contributions. In practice, dispersion relations and 
xPT have yet to yield entirely model-independent predic- 
tions, though considerable insight has been derived from 
their application. 

The benchmark dispersion relation analysis was per- 
formed by Jaffe |3(|| , who employed the narrow resonance 
approximation for the hadronic spectrum that is consis- 
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tent with the large A'";, hmit of QCD and an ansatz for 
the high-Q^ behavior of consistent with known quark 
counting rules. Use of such an ansatz, which introduces a 
degree of model-dependence , was needed to account 
in a physically realistic way for incalculable contributions 
from the higher-mass part of the QCD spectrum. The 
results, which drew upon the measured isoscalar electro- 
magnetic form factors and the flavor content of the lowest 
lying mesons, suggested that G%[ would be both siz- 
able and negative. An updated version of this analysis 
was carried out by the authors of [s^l , who included log- 
arithmic corrections to the power-law asymptotic behav- 
ior used in 'so'l but found no substantial modification of 
the original Jaffe prediction. The results in Eq. how- 
ever, tend to disfavor such a substantially negative value, 
thereby sug gesting a re-examination of the assumptions 
used in [30j. 

Subsequent dispersion theoretic work in [sj |3 Issf 

avoided a priori reliance on the narrow resonance approx- 
imation and high-Q^ ansatz and instead constructed the 
isoscalar and strange magnetic spectral functions from 
measured K-N scattering amplitudes and e"'"e~ partial 
widths. This analysis, which entailed a truncation of the 
spectrum at ^^1 GeV due to the lack of experimental 
input, produced results consistent with the Jaffe predic- 
tions for the low-mass spectral content of G%f and demon- 
strate that the nucleon's "kaon cloud" is dominated by 
the (/)(1020) resonance. Thus, the narrow resonance ap- 
proximation for the low-mass spectral function, as as- 
sumed in [303, appears justified. By itself, however, the 
low-mass spectral content would imply fi^ ~ —0.3, so in- 
clusion of higher-mass contributions to the strangeness 
vector spectral functions appears necessary in order to 
account for the SAMPLE value for G^j. At present, 
no model-independent, first principles treatment of this 
higher-mass region has yet been achieved, though a num- 
ber of one- loop model calculations [sHIs^ls^ls^ have in- 
dicated the importance of the higher-mass states in mov- 
ing the value /is in the right direction to better agree 
with the SAMPLE result. 

More generally, one-loop computations using a 
hadronic basis represent a perturbative approximation 
to the full dispersion theoretic treatment. Such loop cal- 
culations entail truncation of an expansion in the strong 
hadronic coupling g at second order and, therefore, tend 
to be subject to a substantial degree of model-dependent 
ambiguities particular, the one-loop ampli- 

tudes rely on a unitarity violating approximation to pseu- 
doscalar meson-nucleon scattering amplitudes, the omis- 
sion of higher-order (in g) rescattering contributions that 
restore unitarity and generate the physically important 
resonant behavior of the spectral functions, and use of 
unphysical, point-like pseudoscalar vector current form 
factors [4a, |43l ■ In principle, these deficiencies are reme- 
died in a systematic way using xPT, wherein the a pri- 
ori unknown low-energy constants (LEC's) determined 



by fits to experimental data embody the physics omitted 
from the one-loop graphs. In the case of G|^, however, 
this program encounters an intrinsic limitation due to the 
symmetry properties of the strangeness vector current 
(see below), so any one- loop predictions for necessar- 
ily entail unquantifiable model-dependent uncertainties. 

The earliest hadron loop calculations 0> Q> sug- 
gested that the non-resonant part of the nucleon's kaon 
cloud should make a fairly small contribution to G%j, 
an insight confirmed by the subsequent dispersion the- 
ory analyses discussed above. A variation on this theme 
was carried out by Geiger and Isgur [s^, who summed 
up a complete tower of meson-baryon one loop graphs, 
using the quark model to determine the relevant hadronic 
vertices, and found a pattern of cancellations among suc- 
cessively higher-mass intermediate states. Although this 
computation suffers from the same ambiguities as the 
earlier one-loop calculations, it is nevertheless sugges- 
tive that higher mass contributions may, as indicated by 
Eq. 10) and the dispersion theoretic studies, play an im- 
portant role in the dynamics of G\j. 

In general, the systematic, model-independent treat- 
ment of hadronic loop effects in xPT does not yield pre- 
dictions that are independent of the G|^ measurements 
since the operator sj^s contains an SU(3)-singlet com- 
ponent 0- Apart from the SAMPLE result for Gfj 
itself, there exists no other experimental information on 
the SU(3)-singlct component of the nucleon vector cur- 
rent that would allow one to determine the LEC's rel- 
evant to strange magnetism. An exception occurs for 
the strange magnetic radius, < r^ >m, that governs the 
slope of G\j at the origin, for which a parameter-free 
prediction can be made at 0{p^) However, inclu- 

sion of O(p^) loop contributions, nominally suppressed 
by one power of mx/^x where A^^^ = 477,?^ is the chi- 
ral scale, nearly cancels most of the 0{p^) term, leaving 
a residual dependence on an a priori unknown strange 
magnetic radius LEC, 5^ E3|- Rigorously speaking, the 
latter must be determined from future measurements of 
the Q^-dependence of G%[. However, a reasonable range 
may be estimated by comparing the model-independent 
dispersion relation and lattice QCD calculations as in- 
dicated above, leading to — 1 < 5^ < 1. In this case, 
xPT provides reasonable guidance for an extrapolation 
of G%j{Q^) to the photon point. The static moment can 
be written as 

Ais = Gl,(Q2 = o.l)-0.13&^ (6) 

resulting in iig — 0.37 ± 0.20 ± 0.26 ± 0.15, where now 
the two sources of theoretical uncertainty (electroweak 
radiative corrections and strange magnetic radius) have 
been combined. 

Nucleon model calculations, though less transparently 
connected to QCD, have also provided some insights into 
the dynamics of strange quarks. One model that has re- 
ceived considerable attention recently is the chiral quark 
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soliton model (xQSM), in which constituent quarks in- 
teract with the Goldstone bosons of the spontaneously- 
broken, approximate chiral symmetry of QCD in a self- 
consistent way 0, I49II. A recent xQSM computation 
of reported in js^l yields a small positive result, 
0.05<GX^(Q2 ^ 0.1) < 0.1 that is consistent with the 
SAMPLE result. The same model, however, under- 
predicts the proton and neutron magnetic moments by 
~ 40%, and the prediction for is strongly dependent 
on ad hoc assumptions about the long-distance behavior 
of the Goldstone boson field. The computation is nev- 
ertheless suggestive that, at the microscopic level, the 
topology of the QCD vacuum plays a non-trivial role in 
the dynamics of the ss sea. 

Completion of the SAMPLE program and the corre- 
sponding insights derived from a comparison of the ex- 
perimental results to the last decade of theoretical ef- 
forts represents a milestone in the quest to understand 
the quark substructure of the nucleon. Nonetheless, 
work is on-going, and new experimental data will be- 
come available from Jefferson Laboratory and from the 
Mainz PVA4 program 0. An independent determina- 
tion of G%j at low momentum transfer is planned by the 
HAPPEX collaboration using combined forward angle 
parity-violation measurements on hydrogen and helium 
targets 0, |^ . Precise determination of the depen- 
dence of both G% and G|j over the range 0.3 < < 1.0 
(GeV/c)^ will be available from the planned program of 
measurements by the GO collaboration P| , for which data 
taking will soon begin. Theoretically, new efforts to carry 
out unquenched lattice computations and to resolve the 
numerical sampling challenges are underway. One ex- 
pects these efforts to flesh out the framework that has 
now emerged after more than a decade of experimental 
and theoretical work. Of course, future surprises are al- 
ways a possibility, and one may ultimately find that some 
deeper principle governs the dynamics of sea quarks [s^ 
than is apparent to our current understanding. 
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(ANL). 



* Electronic address: beiseQphysics.umd.edu 
[1] D. Kaplan and A. Manohar, Nucl. Phys. B 310, 527 
(1988). 

[2] R.D. McKeown, Phys. Lett. B219, 140 (1989). 

[3] D.H. Beck, Phys. Rev. D39, 3248 (1989). 

[4] B.A. Mueller et al, Phys. Rev. Lett. 78, 3824 (1997). 

[5] D.T. Spayde et al, Phys. Rev. Lett. 84, 1106 (2000). 

[6] R. Hasty et al, Science, 290, 2021 (2000). 

[7] K. Aniol, et al, Phys. Lett B 509, 211 (2001). See also 



K. Aniol, et al, Phys. Rev. Lett. 82, 1096 (1999). 
[8] Mainz experiment 'PVA4", D. von Harrach, spokesman. 
[9] JLAB experiment EOO-006 and EOl-116 ("GO"), 

D.H. Beck, contact. 
[10] JLAB experiment E99-115, K. Kumar and D. Lhuillier, 

contacts. 

[11] S.-L. Zhu, S. J. Puglia, B. R. Holstein, and 

M. J. Ramsey-Musolf, Phys. Rev. D 62, 033008 (2000). 
[12] T.M. Ito, et al, submitted to Phys. Rev. Lett., 

nucl-ex/0310001 
[13] M. J. Musolf, et al, Phys. Rep. 239, 1 (1994). 
[14] See, for example, B.W. Filippone and X. Ji, 

Adv. Nucl. Phys. 26, 1 (2001). 
[15] I. Zel'dovich, JETP Lett. 33, 1531 (1957). 
[16] M. J. Musolf and B. R. Holstein, Phys. Lett B, 242, 461 

(1990). 

[17] GEANT, CERN program library. 

[18] T.V. Kuchto and N.M. Shumeiko, Nucl. Phys. B 219, 
412 (1983). 

[19] L.W. Mo and Y.S. Tsi, Rev. Mod. Phys. 41, 205 (1969). 

[20] J.-W. Chen and X. Ji, Phys. Lett. B 501, 209 (2001). 

[21] J.C. Bergstrom, Phys. Rev. C58,2574 (1998). 

[22] E. Korkmaz, et al, Phys. Rev. Lett. 83, 3609 (1999). 

[23] See also D.T. Spayde, Ph.D. thesis, Univ. Mary- 
land, 2001 (unpublished but available at 
http: / /www. physics. umd.edu /enp /theses / I . 

[24] S. Galster, et al, Nucl. Phys. B 32, 221 (1971). 

[25] D.H. Beck and B.R. Holstein, Int. J. Mod. Phys. E 10, 
1 (2001). 

[26] S. J. Dong and K. F. Liu, Phys. Lett. B 328, 130 (1994). 
[27] S.J. Dong, K.F. Liu and A.G. Wilhams, Phys. Rev. D 

58, 074504, (1998). 
[28] R. Lewis, W. Wilcox, and R.M. Woloshyn, Phys. Rev. D 

67, 013003 (2003). 
[29] D. Leinweber and A.W. Thomas, Phys. Rev. D 62, 

074505 (2000). 
[30] R.L. Jaffe, Phys. Lett. B229, 275 (1989). 
[31] H. Forkel, Phys. Rev. C 56, 510 (1997). 
[32] H. W. Hammer, U. G. Meissner and D. Drechsel, Phys. 

Lett. B 367, 323 (1996). 
[33] M.J. Ramsey-Musolf and H.-W. Hammer, 

Phys. Rev. Lett. 80, 2539 (1998). 
[34] H.-W. Hammer and M.J. Ramsey-Musolf, 

Phys. Rev. C60, 045204 (1999) [Erratum-ibid. C 

62, 049902 (2000)]. 
[35] H.-W. Hammer and M.J. Ramsey-Musolf, 

Phys. Rev. C60, 045205 (1999) [Erratum-ibid. C 

62, 049903 (2000)]. 
[36] H. W. Hammer and M. J. Ramsey-Musolf, Phys. Lett. B 

416, 5 (1998). 

[37] L. L. Barz, H. Forkel, H. W. Hammer, F. S. Navarra, 
M. Nielsen and M. J. Ramsey-Musolf, Nucl. Phys. A 640, 
259 (1998). 

[38] P. Geiger and N. Isgur, Phys. Rev. C55, 299 (1997). 
[39] D. O. Riska, Few Body Syst. Suppl. 10, 415 (1999). 
[40] H. Ito and M. J. Ramsey-Musolf, Phys. Rev. C 58, 2595 
(1998). 

[41] T.R. Hemmert, U-G. Meissner, and S. Stieninger, 

Phys. Lett. B437, 184 (1998). 
[42] H.-W. Hammer, S.J. Puglia, M.J. Ramsey-Musolf, and 

S.-L. Zhu, Phys. Lett. B 562, 208 (2003). 
[43] W. Koepf, E. M. Henley and S. J. Pollock, Phys. Lett. B 

288, 11 (1992). 

[44] M. J. Musolf and M. Burkardt, Z. Phys. C 61, 433 (1994). 



6 



[45] H. Forkel, M. Nielsen, X. m. Jin and T. D. Cohen, Phys. 

Rev. C 50, 3108 (1994). 
[46] P. Federbush, M.L. Goldberger, and S.B. Treiman, Phys. 

Rev. 112, 642 (1958). 
[47] M. J. Musolf, H. W. Hammer and D. Drechsel, Phys. Rev. 

D 55, 2741 (1997) [Erratum-ibid. D 62, 079901 (2000)]. 
[48] C.V. Christov, et al, Prog. Part. Nucl. Phys. 37, 91 

(1996). 



[49] R. Alfoker, H. Reinhardt, and H. Weigel, Phys. Rep. 265, 
139 (1996). 

[50] A. Silva, H-C Kim, and K. Goeke, 

|arXiv:h ep-ph/0210189 
[51] JLab experiment EOO-114, D. Armstrong and 

R. Michaels, contacts. 
[52] X. Ji and J. Tang, Phys. Lett. B 362, 182 (1995). 



